bep-lat/02030T3 



The topological susceptibility of 
SU(3) gauge theory near T c 

Christof Gattringer^, Roland Hoffmann and Stefan Schaefer 

Institut fur Theoretische Physik, Universitat Regensburg 
D-93040 Regensburg, Germany 



Abstract 

We compute the topological susceptibility Xt in SU(3) lattice gauge 
theory using fermionic methods based on the Atiyah-Singer index theo- 
rem. Near the phase transition we find a smooth crossover behavior for 
Xt with values decreasing from (191(5) MeV) 4 to (100(5) MeV) 4 as we 
increase the temperature from 0.88 T c to 1.31 T c , showing that topo- 
logical excitations exist far above T c . Our study is the first large scale 
analysis of the topological susceptibility at high temperature based 
on the index theorem and the results agree well with field theoretical 
methods. 
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Relating topology to the dynamics of a gauge theory is a fascinating 
idea which can provide deep insights to the non-perturbative structure of 
the theory. For QCD an object of particular interest is the topological 
susceptibility xt defined as 

xt = ^(Q 2 )- (i) 

Q denotes the topological charge obtained e.g. as an integral over F^ U F^ U . 
For QCD xt plays an important role in the low energy behavior where e.g. the 
so-called Witten-Veneziano formula [j^] relates the quenched xt to the meson 
masses. For pure gauge theory xt is a measure for the net abundance of 
topological fluctuations and can provide insight to instanton liquid models 
||]. Here the behavior around the phase transition is of particular interest, 
since new machines such as RHIC are starting to probe hot QCD. xt is 
expected to show a smooth crossover behavior into the high temperature 
phase as topological excitations die out. 

The computation of xt near T c is an inherently non-perturbative calcula- 
tion and can e.g. be attacked in lattice gauge theory. While for SU(3) lattice 
gauge theory at zero temperature there are several recent calculations 
only a single modern study || of the critical behavior near T c can be found 
in the literature (for an early attempt see ||). In this article we improve 
on this situation and present a systematic study of xt near the critical tem- 
perature with high statistics and an assessment of finite size effects from a 
comparison of data from different lattice sizes. 

The traditional approach || to the calculation of the topological charge 
was a discretization of F^ U F^ V (the so-called field theoretical method). Lat- 
tice discretizations of F^F^ are very sensitive to ultraviolet fluctuations 
and can be used only after some cooling or smearing procedure has been 
applied to the gauge field. Such a cooling step, however, has the potential 
to destroy topological lumps and to alter the outcome for xt- Thus it is 
important to check whether there is consistency with alternative methods. 

In this letter we compute the topological charge using a fermionic method 
(compare also [f|, 0). For the continuum Dirac operator D the Atiyah- 
Singer index theorem |@] relates the topological charge Q to the index of D 
which can be written as the difference of the numbers of left-handed and 
right-handed zero- modes. For the standard Wilson-Dirac operator such a 
determination of Q through the zero-modes fails for numerically accessible 
values of the cutoff since the would-be zero-modes mix with the correspond- 
ing modes in the doubler branches ||. 



1 



With the re-discovery of the Ginsparg- Wilson relation [[!(]] for lattice 
Dirac operators it was understood how to formulate chirally symmetric 
fermions on the lattice. In particular Ginsparg- Wilson fermions give rise 
to an index theorem [11] already at finite cutoff. It was found 0, HI that 

Q = n _ - n+ = ^Tt l5 D = ^ J (PxF^F^x) + 0(a 2 ) . 

(2) 

In the first line of the equation n+ (n_) denotes the number of zero- modes 
with positive (negative) chirality. The problem with the doublers is resolved 
since their would-be zero-eigenvalues are shifted to 2/ a. The great disad- 
vantage of an exact solution of the Ginsparg- Wilson equation such as the 
overlap operator [ 13 1 is the high cost of a numerical implementation. 

A possible, considerably cheaper approach is to use an approximate so- 
lution of the Ginsparg Wilson equation such as a finite parameterization of 
the fixed point action [14]. In this study we use chirally improved fermions 
[15 1 which arise from a systematic expansion of a solution of the Ginsparg- 
Wilson equation. They are sufficiently chiral to allow for an unambiguous 
identification of the zero-modes but are numerically considerably cheaper 
than overlap fermions and allow to obtain the statistics necessary for a pre- 
cise measurement of xt- 

Our (quenched) SU(3) gauge configurations were generated with the 
Liischer-Weisz action jl6|. The leading coupling 0i varied in a range from 
Pi = 8.10 to Pi = 8.60. When using the Sommer parameter ro = 0.5 fm 
to set the scale this corresponds [17| to a range of lattice spacings from 
a = 0.084 fm to a = 0.125 fm (compare Table [l] below). For details con- 
cerning the Monte Carlo see ||18| . We used two settings of lattices: Firstly, 
L 4 lattices with L = 16 and L = 12. These ensembles are all in the low 
temperature phase of QCD, where quarks are confined and chiral symmetry 
is broken. Secondly, we use lattices of size 6 x L 3 with L = 20, L = 16 
and L = 12. These latter ensembles give rise to temperatures ranging from 
values below the critical temperature (0.88 T c ) to values above the phase 
transition (1.31 T c ). The details for the parameters of our runs can be found 
in Table |l[ For the Liischer-Weisz action the critical temperature was com- 
puted in []l9| using Polyakov loops and also the gap of the spectrum of the 
Dirac operator. Both the chiral transition and the deconfinement transition 
were found to have a critical temperature of T c = 300(3) MeV. 

Since the chirally improved fermions we are using here are only an ap- 
proximation of a solution of the Ginsparg- Wilson equation a few comments 
on our determination of the topological charge are in order here. One of the 
effects of using an approximation of a Ginsparg- Wilson operator are fluctu- 
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lattice 


HI 


a [fm] 


T [MeV] 


N t 

x Y con j s 


i 

yt [Me VI 

At L v J 


6 x 20 3 


8.10 


0.125(1) 


264(2) 


400 


194(5) 


6 x 20 3 


8.20 


0.115(1) 


287(3) 


400 


185(5) 


6 x 20 3 


8.25 


0.110(1) 


299(3) 


400 


162(5) 


6 x 20 3 


8.30 


0.106(1) 


311(3) 


400 


144(5) 


6 x 20 3 


8.45 


0.094(1) 


350(4) 


400 


124(5) 


6 x 20 3 


8.60 


0.084(1) 


392(5) 


400 


93(5) 


6 x 16 3 


8.10 


0.125(1) 


264(2) 


800 


200(4) 


6 x 16 3 


8.20 


0.115(1) 


287(3) 


800 


186(4) 


6 x 16 3 


8.25 


0.110(1) 


299(3) 


800 


170(4) 


6 x 16 3 


8.30 


0.106(1) 


311(3) 


800 


153(4) 


6 x 16 3 


8.45 


0.094(1) 


350(4) 


800 


123(4) 


6 x 16 3 


8.60 


0.084(1) 


392(5) 


800 


105(5) 


6 x 12 3 


8.10 


0.125(1) 


264(2) 


1200 


191(4) 


6 x 12 3 


8.20 


0.115(1) 


287(3) 


1200 


183(4) 


6 x 12 3 


8.25 


0.110(1) 


299(3) 


1200 


170(4) 


6 x 12 3 


8.30 


0.106(1) 


311(3) 


1200 


153(4) 


6 x 12 3 


8.45 


0.094(1) 


350(4) 


1200 


118(4) 


6 x 12 3 


8.60 


0.084(1) 


392(5) 


1200 


101(5) 


16 4 


8.10 


0.125(1) 


99(1) 


200 


185(6) 


16 4 


8.20 


0.115(1) 


107(1) 


200 


194(7) 


16 4 


8.30 


0.106(1) 


117(1) 


200 


189(6) 


16 4 


8.45 


0.094(1) 


131(1) 


200 


194(6) 


16 4 


8.60 


0.084(1) 


147(2) 


200 


192(7) 


12 4 


8.10 


0.125(1) 


132(1) 


400 


200(5) 


12 4 


8.30 


0.106(1) 


155(1) 


400 


188(5) 


12 4 


8.45 


0.094(1) 


175(2) 


400 


191(5) 



Table 1: Parameters and results. We give the size of our lattices, the in- 
verse gauge coupling f3±, the corresponding lattice spacing a in fermi, the 
temperature in MeV, the number N con f s of configurations in our ensembles 
and our results for the topological susceptibility. 
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Figure 1: Distribution of the real eigenmodes of the chirally improved Dirac 
operator for the (3 = 8.20 ensembles. The dashed histogram represents the 
data from the 16 4 lattice while the full line is for 6 x 20 3 . 



ations of the zero-eigenvalues. Instead of being located exactly at the origin 
the zero-eigenvalues manifest themselves as small real eigenvalues. It can be 
demonstrated [2C] that for discretized instantons the size of these would-be 
zero-eigenvalues increases as the radius of the instanton decreases. Simi- 
larly, also for the thermalized configurations one finds a distribution of real 
eigenvalues. For two of our ensembles (16 4 and 6 x 20 3 both at fa = 8.20) 
we show this distribution in Fig. |]. The distribution of the position x of 
the real eigenvalues has a pronounced peak at the origin and a tail which 
extends towards larger values. The tail comes from very small instantons, 
so-called dislocations. 

We compute our eigenvalues using the implicitly restarted Arnoldi method 
plf . This algorithm computes a given number N of eigenvalues ordered with 
respect to their absolute value. We always compute the N = 50 smallest 
eigenvalues, except for the 16 4 lattices where we have ./V = 30. A potential 
source of error is a too low value of N. Sufficiently many eigenvalues have 
to be computed in order to capture all of the peak in the distribution of 
Fig. |l[ As is obvious from Fig. [l] with our choice of N we obtain all of the 
peak and a large portion of the tail. Setting a cut for the eigenvalues at 
e.g. 0.1 leaves the results for \t unchanged showing that dislocations make 
up only a small contribution to the topological susceptibility. We remark 
that there is no danger of mistaking a would-be zero eigenvalue for an eigen- 
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value with a small but non- vanishing imaginary part. A few lines of algebra 
p2f show that for a 75-hermitian Dirac operator the matrix element V^75^ 
of an eigenvector ip is non- zero only for eigenvectors with real eigenvalues. 
For our approximation of a Ginsparg- Wilson Dirac operator we find values 
of 1^ TsV'l ranging from 0.8-0.9 in the peak of the distribution down to 0.4 
for the very end of the tail. 

We summarize our procedure for determining Q as follows: We compute 
the 50 (30 for 16 4 ) eigenvalues of the Dirac operator and for the eigenvectors 
with real eigenvalues evaluate ip^j^ijj. We take n + (n_) to be the number 
of modes where V^75^ is positive (negative). The topological charge is then 
computed as in Eq. @) by Q = n„ — n + . 

For a subsample of 200 configurations on the 12 lattices at Pi = 8.10 and 
j3\ = 8.45 we cross-checked our procedure with the results obtained from the 
overlap operator. For the overlap operator the zero-modes have eigenvalues 
exactly at zero since also modes from dislocations are projected onto exact 
zero-modes. The corresponding 75 matrix elements are always ip'^ip = ±1. 
When comparing the determination of Q for the individual configurations 
we found a discrepancy for 2% for the configurations at f3\ = 8.45 and 9% 
for (5\ = 8.10. This demonstrates that as one goes closer to the continuum 
limit the two definitions agree better. Furthermore, the difference in Q was 
always only one unit such that the results for xt agree surprisingly well: For 
the subsamples of 200 configurations we found at /?i = 8.45 a value of xt = 
(182 (7) MeV) 4 for both methods. At = 8.10 we obtained a value of xt = 
(196 (7) MeV) 4 with the overlap operator and xt = (197 (7) MeV) 4 with the 
chirally improved operator. Our test shows that the two methods agree very 
well, however, using the overlap operator is by a factor of 10 more costly 
than working with the chirally improved operator. We remark that a direct 
comparison of the field theoretical approach and the fermionic method on 
single configurations can be found in 

We begin with a discussion of the results for \t as a function of j3\. We 
plot the corresponding data in Fig. ^. For time extent 6 the critical value of 
f}\ was determined [19] to be (3\ = 8.26(2). We mark the transition with a 
dashed vertical line in Fig. ||. For the 6 x L 3 configurations we used values 
of 0i ranging from 8.10 to 8.60 giving rise to ensembles on both sides of 
the transition. The corresponding data are represented in Fig. || by squares 
(6 x 12 3 ), triangles (6 x 16 3 ) and diamonds (6 x 20 3 ). We connect the 
symbols to guide the eye. When comparing the results for different volumes 
we find that the data agree well within error bars. This shows that for the 
chosen lattice sizes we do not encounter finite size problems. For /?i = 8.10 
our values for xt scatter around (190 MeV) 4 , which is the value in the low 
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Figure 2: The topological susceptibility as a function of ft We show our 
results for the zero temperature ensembles on 12 4 (circles) and 16 4 (filled 
circles) as well as for the finite temperature ensembles on 6 x 12 3 ,6 x 16 3 
and 6 x 20 3 lattices (squares, triangles, diamonds). The dashed vertical line 
marks the phase transition. 

temperature phase. As one increases ft, xt starts to drop already before 
the critical value ft = 8.26 is reached. The slope of the curve is largest near 
the critical ft. The decrease slows down as ft is increased further and xt 
reaches a value of xt ~ (100 MeV) 4 at ft = 8.60. 

The decreasing curve for xt from the ensembles at high temperature can 
now be compared to the low temperature data with the same values of ft, 
i.e. the same cutoff. The corresponding data are represented by open circles 
(12 4 ) and filled circles (16 4 ). Again we find that the data from the two 
different volumes are well consistent within error bars and we do not face 
finite size problems. The numbers for xt remain near (190 MeV) 4 for all 
values of ft. 

It is interesting to combine the results from the zero temperature en- 
sembles and compare them to the results available in the literature. The 
combined result from our 16 4 lattices is xt = (191(5) MeV) 4 . The more 
recent results in the literature range between xt = (175(5) MeV) 4 to xt = 
(203(5) MeV) 4 . This demonstrates that our zero temperature result for the 
Liischer-Weisz action determined from the index theorem agrees well with 
published data. 

Let us now return to the more interesting behavior near the phase tran- 
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Figure 3: The topological susceptibility as a function of the temperature. 
We combine our results from all lattices: 12 4 , 16 4 , 6 x 12 3 , 6 x 16 3 and 6 x 20 3 
(circles, filled circles, squares, triangles, diamonds). The dashed vertical line 
marks the phase transition. 

sition. In Fig. || we present our data as a function of the temperature. We 
include also our data for the L 4 lattices in order to have the baseline in 
the low temperature phase. The data from our high temperature, 6 x L 3 
configurations cover a temperature range from 264 MeV to 392 MeV. Again 
we mark the phase transition by a vertical dashed line. As already men- 
tioned above, the data from different volumes agree very well and we do 
not encounter finite size problems. The topological susceptibility starts to 
deviate from its zero temperature value (xt ~ (190 MeV) 4 ) at a temperature 
of T ~ 275 MeV which is still below the deconfinement transition. The drop 
continues and becomes steepset at T c and then leans back slightly. At our 
largest temperature T = 392 MeV (= 1.31 T c ) the decrease is still substan- 
tial and xt has reached xt ~ (100 MeV) 4 which is about half of its zero 
temperature value. A linear extrapolation of the decrease of xl^ gives a 
temperature of about 600 MeV where xt becomes compatible to zero. This 
would be about twice the critical temperature, but since the decrease of 
x]^ 4 is certainly not linear throughout, \t ~ is presumably reached at a 
temperature even higher than 600 MeV. Our data are in reasonably well 
agreement with the other result for xt near T c available in the literature || . 

In this letter we have performed the first large scale study of the finite 
temperature behavior of xt based on the index theorem. Fermionic methods 
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provide an alternative approach free from possible ambiguities due to cool- 
ing or smearing. From a comparison of data on different volumes we find 
that finite size effects are under good control. We find that the topological 
susceptibility extends into the high temperature phase and some topological 
excitations can be found up to 2 T c . 
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